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ABSTRACT 


A  set  of  multiscale,  nested,  idealized  numerical  simulations  of  mesoscale  convective  systems  (MCSs)  and 
mesoscale  convective  vortices  (MCVs)  was  conducted.  The  purpose  of  these  simulations  was  to  investigate 
the  dependence  of  MCV  development  and  evolution  on  background  conditions  and  to  explore  the  relation¬ 
ship  between  MCVs  and  larger,  moist  baroclinic  cyclones.  In  all  experiments,  no  mesoscale  convective  system 
(MCS)  developed  until  a  larger-scale,  moist  baroclinic  system  with  surface  pressure  amplitude  of  at  least 
2  hPa  was  present.  The  convective  system  then  enhanced  the  development  of  the  moist  baroclinic  system  by 
its  diabatic  production  of  eddy  available  potential  energy  (APE),  which  led  to  the  enhanced  baroclinic 
conversion  of  basic-state  APE  to  eddy  APE.  The  most  rapid  potential  vorticity  (PV)  development  occurred 
in  and  just  behind  the  leading  convective  line.  The  entire  system  grew  upscale  with  time  as  the  newly  created 
PV  rotated  cyclonically  around  a  common  center  as  the  leading  convective  line  continued  to  expand  outward. 
Ten  hours  after  the  initiation  of  deep  moist  convection,  the  simulated  MCV  radii,  heights  of  maximum  winds, 
tangential  velocity,  and  shear  corresponded  reasonably  well  to  their  counterparts  in  BAMEX.  The  increasing 
strength  of  the  simulated  MCVs  with  respect  to  larger  values  of  background  CAPE  and  shear  supports  the 
hypothesis  that  as  long  as  convection  is  present,  CAPE  and  shear  both  add  to  the  strength  of  the  MCV. 


1.  Introduction 


PV  concepts  to  explore  the  relationship  between  the 
balanced  circulation  of  the  vortex  and  its  interaction 
with  the  vertical  shear  in  which  the  vortex  resides.  The 
same  conceptual  model  has  been  invoked  (Fritsch  et  al. 
1994;  Galarneau  et  al.  2009)  to  explain  the  persistence, 
cyclic  redevelopment,  and  upscale  growth  of  MCVs.  In 
this  model,  vertical  gradients  of  diabatic  heating  pro¬ 
duce  a  PV  maximum  in  the  lower  to  mid  troposphere 
overlain  by  a  PV  relative  minimum  near  the  tropopause. 
Although  the  PV  structure  resides  in  a  layer  of  vertical 
shear  (due  to  thermal  wind  balance)  that  would  tend  to 
destroy  the  vortex  over  time,  the  balanced  cyclonic  mo¬ 
tion  associated  with  the  vortex  brings  about  a  region  of 
warm,  moist  advection  and  lift  on  the  downshear  side  of 
the  vortex.  The  lift  brings  about  further  convection  and 
diabatic  heating  that  helps  maintain  the  vortex  in  the 
face  of  the  sheared  environment. 


The  development  of  mesoscale  convective  vortices 
(MCVs)  is  one  example  of  a  process  in  which  diabatic 
heating  plays  a  central  role  in  the  intensification  and 
maintenance  of  an  atmospheric  vortex  (Fritsch  et  al.  1994; 
Raymond  and  Jiang  1990;  Jiang  and  Raymond  1995; 
Houze  2004;  Hertenstein  and  Schubert  1991;  Davis  and 
Weisman  1994).  Using  potential  vorticity  (PV)  concepts, 
the  intensification  can  be  explained  by  the  vertical  gra¬ 
dient  of  diabatic  heating  within  the  convective  system. 
Raymond  and  Jiang  (1990)  and  Jiang  and  Raymond 
(1995)  developed  a  nonlinear  balance  model  using  these 
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The  same  PV  concepts  have  been  employed  as  a  pro¬ 
posed  mechanism  for  the  initial  development  of  moist 
baroclinic  systems  (Montgomery  and  Farrell  1991;  Moore 
and  Montgomery  2004,  2005,  hereafter  MM04  and 
MM05).  In  the  moist  baroclinic  case,  the  diabatic  heat¬ 
ing  commences  in  a  region  of  lift  well  downshear  from 
the  initiating  midlevel  disturbance,  and  the  initial  vortex 
that  develops  is  resistant  to  the  effects  of  the  shear  be¬ 
cause  of  the  reinforcing  effects  of  the  latent  heat  release 
in  clouds  on  its  downshear  side.  MM05  coined  the  term 
“diabatic  Rossby  vortex”  for  such  vortical  disturbances, 
although  in  an  earlier  study  Snyder  and  Lindzen  (1991) 
referred  to  the  wavelike  counterparts  as  diabatic  Rossby 
waves. 

Whereas  cold  season  baroclinic  cyclones  ultimately 
derive  a  majority  of  their  energy  from  the  available 
potential  energy  (APE)  of  the  basic  state  (Davis  et  al. 
1993),  the  diabatic  production  of  PV  can  be  a  strong 
contributor  to  their  development,  especially  over  the 
oceans  (Montgomery  and  Farrell  1992;  Whitaker  and 
Davis  1994;  MM04;  MM05;  Moore  et  al.  2008).  In  warm 
season  situations,  the  background  APE  may  be  some¬ 
what  less,  but  the  convective  instability  also  allows  locally 
stronger  diabatic  heating;  like  their  cold  season  counter¬ 
parts,  the  warm  season  disturbances  move  along  larger- 
scale  temperature  gradients.  Additionally,  MM04  show 
a  region  of  scale  independence  of  moist  baroclinic  growth 
below  a  certain  horizontal  wavelength,  motivating  an 
extension  of  these  studies  into  environments,  more  typi¬ 
cal  of  warm  season  conditions  (Conzemius  et  al.  2007, 
hereafter  Cl)  that  allow  the  development  of  some  CAPE. 

In  the  idealized  simulations  of  Cl,  it  was  found  that 
the  development  of  organized,  deep,  moist  convection 
required  the  presence  of  larger-scale  forcing.  The  de¬ 
velopment  of  the  mesoscale  convective  system  (MCS) 
occurred  in  conjunction  with  the  deepening  of  the  most 
rapidly  growing  moist  baroclinic  wave,  regardless  of  any 
attempt  to  minimize  the  effects  of  the  larger  moist  baro¬ 
clinic  system  in  the  simulations.  The  overall  conclusion 
was  that  an  important  scale  interaction  existed,  in  which 
the  development  of  the  MCS  was  dependent  on  the 
deep,  persistent  upward  vertical  motion  and  the  moist, 
low-level  inflow  provided  by  the  larger-scale  system. 

In  the  present  study,  another  avenue  we  wish  to  ex¬ 
plore  is  the  reverse  interaction — whether  the  MCVs  have 
a  direct  feedback  on  the  growth  rate  of  synoptic-scale, 
baroclinic  waves.  Any  vortex  that  has  a  localized  distri¬ 
bution  in  physical  space  will  project  onto  a  spectrum  of 
scales  that  will  include  large-scale  baroclinic  waves  (Cl), 
which  may  or  may  not  be  exponentially  unstable.  An 
MCV  developing  within  a  moist  baroclinic  system  should 
be  no  exception.  Additionally,  the  energy  conversion 
calculations  in  Cl  revealed  that  although  the  diabatic 


production  of  eddy  APE  is  a  significant  contributor  to 
the  net  eddy  APE,  the  conversion  of  mean  state  APE  to 
eddy  APE  (baroclinic  conversion)  increased  in  tandem 
with  the  diabatic  production  of  APE.  This  result  sug¬ 
gests  an  important  upscale  feedback  from  the  mesoscale 
(MCV/MCS)  to  synoptic  scales  in  the  development  of 
moist  baroclinic  cyclones.  Although  less  baroclinically 
unstable  environments  can  be  expected  in  the  summer 
when  MCVs  are  much  more  prevalent,  it  can  still  be 
expected  that  the  MCVs  play  a  role  in  the  maintenance 
of  synoptic-scale  systems  in  the  midlatitudes  in  summer, 
particularly  since  they  often  appear  to  undergo  an  upscale 
growth  (Fritsch  et  al.  1994;  Cl;  Galarneau  et  al.  2009).  In 
more  weakly  sheared  environments,  there  may  also  be 
some  implications  for  the  development  and  maintenance 
of  easterly  waves  and  the  mesoscale  convective  systems 
embedded  within  them  as  recently  hypothesized  in  the 
marsupial  theory  of  tropical  cyclogenesis  (Dunkerton 
et  al.  2008;  see  their  hypothesis  3). 

A  second  goal  of  the  present  study  is  to  take  a  more 
detailed  look  at  the  PV  development  process  and,  in 
particular,  the  process  by  which  the  MCV  grows  upscale 
to  the  point  that  it  becomes  not  easily  distinguished  from 
the  developing  moist  baroclinic  cyclone  (Galarneau 
et  al.  2009).  One  question  that  arises  is  whether  or  not 
the  upscale  growth  occurs  via  diabatic  production  or  via 
axisymmetrization  of  the  PV. 

In  the  spring  and  summer  of  2003,  the  Bow  Echo  and 
Mesoscale  Convective  Vortex  Experiment  (BAMEX; 
Davis  et  al.  2004)  took  place  with  an  aim  to  observe,  on 
multiple  scales,  the  process  of  MCV  formation  within 
MCSs.  Although  MCVs  have  historically  been  assumed 
to  form  primarily  within  environments  where  the  shear 
is  relatively  small  and  dynamically  of  second  order 
(Bartels  and  Maddox  1991;  Cotton  and  Anthes  1989; 
Houze  2004),  other  studies  have  found  them  to  form  in 
both  weakly  and  strongly  sheared  background  environ¬ 
ments  (Trier  et  al.  2000;  Davis  et  al.  2002).  The  more 
recent  BAMEX  data  demonstrated  MCVs  to  be  present 
in  both  strong  and  weak  deep-layer  shear  environments 
(Davis  and  Trier  2007).  Given  these  results,  a  more  com¬ 
plete  understanding  of  the  MCV  problem  with  varying 
vertical  shear  seems  highly  warranted.  Since  the  Cl  study 
only  covered  background  conditions  of  relatively  modest 
shear,  we  seek  here  to  extend  the  investigation  using 
idealized  simulations  into  more  strongly  sheared  envi¬ 
ronments.  For  similar  reasons,  the  investigations  need  to 
explore  environments  with  larger  CAPE.  The  recent 
results  of  Davis  and  Trier  (2007),  which  examined  MCV 
structure  within  environments  of  varying  CAPE  and 
shear,  also  naturally  motivate  us  to  determine  whether 
similar  structures  can  be  simulated  with  similar  ante¬ 
cedent  conditions  as  in  BAMEX. 
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2.  Methods 

The  idealized  simulations  were  constructed  using  the 
same  basic  methods  as  in  MM05  and  Cl.  A  meridionally 
varying,  Eady-type  basic  state  was  constructed,  and  the 
initial  state  of  the  simulation  was  created  by  inserting  a 
low-level,  warm-core  vortex  with  zero  interior  PV  anom¬ 
alies  into  the  basic  state.  The  fifth-generation  Pennsyl¬ 
vania  State  University-National  Center  for  Atmospheric 
Research  (NCAR)  Mesoscale  Model  (MM5)  was  used  to 
numerically  integrate  the  initial  state  forward  in  time. 

a.  Basic-state  description 

In  Cl,  the  construction  of  the  basic  state  started  by 
specifying  a  zero  surface  pressure  gradient  and  a  surface 
meridional  temperature  gradient  whose  magnitude  varied 
according  to  a  Gaussian  function  of  latitude  centered  on 
45°N.  This  meridional  temperature  gradient  was  held  con¬ 
stant  with  height  in  the  troposphere,  and  thermal  wind 
balance  was  used  to  construct  the  basic  state  zonal  wind. 
These  same  methods  are  used  in  the  present  study. 

Because  the  temperature  lapse  rate  in  a  moist  neutral 
vertical  profile  varies  strongly  as  a  function  of  temper¬ 
ature  (because  of  the  exponential  dependence  of  satu¬ 
ration  vapor  pressure  on  temperature),  Cl  had  to  use 
a  latitude-dependent  LCL  height  to  produce  the  Eady- 
type  shear  in  the  basic  state.  However,  in  the  present 
study,  the  intent  was  to  create  a  background  environ¬ 
ment  of  larger  CAPE,  so  a  constant  temperature  lapse 
rate  of  8  K  km^'  was  used  above  the  PBL  top  through¬ 
out  the  tropospheric  portion  of  the  domain.  To  avoid 
convection  initiating  everywhere  in  the  CAPE-laden 
portion  of  the  domain,  a  temperature  jump  of  1.5  K  and 
a  dewpoint  depression  of  1.5  K  were  used  at  the  top  of 
the  PBL.  The  PBL  top  was  set  to  a  constant  value  of 
850  hPa  in  the  entire  domain.  The  potential  tempera¬ 
ture  at  the  PBL  top  was  the  same  as  the  surface  potential 
temperature  (i.e.,  a  well-mixed  PBL),  and  the  mixing 
ratio  was  also  held  constant  in  the  PBL.  The  PBL  mixing 
ratio  was  calculated  using  the  temperature  and  the  dew¬ 
point  depression  at  the  PBL  top.  An  upper  bound  of 
16  g  kg^^  was  used  to  limit  the  CAPE  in  the  basic  state; 
with  this  limit,  the  magnitude  of  the  PBL  top  dewpoint 
depression  in  the  southern  portion  of  the  domain  was 
larger  than  1.5  K. 

The  tropopause  pressure  was  allowed  to  vary  merid¬ 
ionally  in  the  shape  of  half  the  period  of  a  sine  function 
over  the  width  of  the  baroclinic  zone  with  the  maximum 
slope  at  the  center  of  the  baroclinic  zone  (45°N).  The 
total  width  of  the  baroclinic  zone  was  60  degrees  latitude 
for  most  cases.  The  vertical  profile  of  relative  humidity 
was  that  of  a  hyperbolic  tangent  function  with  the  hu¬ 
midity  decreasing  from  its  value  at  the  top  of  the  PBL 


to  25%  at  the  tropopause,  and  the  stratospheric  relative 
humidity  was  held  constant  at  25%.  A  sample  sound¬ 
ing  from  the  center  of  the  baroclinic  zone  is  shown  in 
Fig.  1. 

b.  Experiments 

We  adopt  a  nomenclature  for  these  simulations  in  the 
format  CAccSHssBZww,  where  “cc”  is  the  surface-based 
CAPE  (through  the  rest  of  the  document,  we  refer  to 
CAPE  as  surface-based  CAPE)  in  the  center  of  the  baro¬ 
clinic  zone  (hundreds  of  J  kg^'^),  “ss”  the  shear  in  the  basic 
state  (tenths  of  m  s^'^  km^'^)  and  “ww”  the  width  of  the 
baroclinic  zone  (degrees  latitude).  The  format  is  designed 
to  allow  the  reader  to  quickly  identify  the  basic-state  con¬ 
figuration  during  the  results  discussion. 

Five  separate  basic-state  configurations  have  been 
constructed  for  the  exploration  of  the  parameter  space 
(see  Table  1).  The  first  is  the  configuration  from  Cl  and 
is  referred  to  as  the  CA00SH15BZ60  simulation.  This 
simulation  has  a  background  profile  that  is  neutral  to  the 
ascent  of  saturated  parcels  (i.e.,  moist  adiabatic)  and  is 
the  starting  point  for  our  exploration  of  the  parameter 
space.  The  second  basic  state  (CA22SH15BZ60)  con¬ 
tains  a  steeper  midlevel  temperature  lapse  rate  to  allow 
larger  CAPE  in  the  initial  state.  It  has  2200  J  kg^^ 
CAPE  in  the  center  of  the  baroclinic  zone  and  the  same 
shear  as  in  the  CA00SH15BZ60  experiment.  To  further 
explore  the  range  of  shear  that  has  been  observed  in 
MCVs  (Trier  et  al.  2000;  Davis  and  Trier  2007),  the  third 
basic  state  (CA15SH30BZ60)  has  similarly  large  CAPE 
but  twice  the  shear  of  the  first  two  simulations.  To  keep 
the  surface  temperature  in  the  southern  portion  of  the 
domain  within  reasonable  atmospheric  limits,  the  cen¬ 
tral  baroclinic  zone  temperature  in  the  CA15SH30BZ60 
simulations  was  set  to  300  K,  reducing  the  CAPE  there 
to  1553  J  kg^'^.  Motivated  by  the  small  values  of  shear 
observed  in  BAMEX  intensive  observing  period  (lOP) 
8  (Davis  and  Trier  2007),  we  also  performed  a  simula¬ 
tion  with  half  the  shear  of  the  CA00SH15BZ60  sim¬ 
ulation  but  large  CAPE  (CA22SH08BZ60).  Finally,  a 
simulation  with  a  narrower  baroclinic  zone  was  con¬ 
ducted  (CA15SH30BZ30),  the  purpose  of  which  was  to 
mimic  the  more  typical  structure  of  warm  season  baro¬ 
clinic  zones,  which  are  often  latitudinally  much  narrower 
than  their  cold  season  counterparts.  The  narrower  baro¬ 
clinic  zone  would  allow  the  MCV  to  move  northward  out 
of  this  zone  more  easily  and  thereby  experience  changes 
in  forcing  with  regard  to  the  thermal  gradient  and  vertical 
shear.  Additionally,  we  wished  to  reduce  the  CAPE  to 
the  north  of  the  baroclinic  zone  since  it  appeared  that  in 
the  other  experiments,  the  CAPE  there  caused  more 
widespread  deep  convective  development  than  is  typi¬ 
cally  found  in  observed  systems. 
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Fig.  1.  Skew-T-log-/?  diagram  in  the  center  of  the  baroclinic  zone  (45°N)  for  the  basic  state  of 
the  CA22SH15BZ60  simulation  (see  Table  1). 


c.  Initial  vortex 

During  the  initial  setup  of  the  CA22SH15BZ60  sim¬ 
ulation,  the  shape  of  the  initial  low-level  vortex  tem¬ 
perature  perturbation  was  set  to  the  second  derivative 
of  the  Gaussian  profile  (the  isolated  vortex),  with  an 
e-folding  radius  of  300  km.  The  reason  for  doing  so  was 
to  minimize  the  initial  projection  on  the  most  rapidly 
growing  and  most  baroclinic  wavenumber.  However, 
since  the  isolated  initial  radial  profile  only  served  to 
delay  the  MCS  development  relative  to  the  simulations 
initialized  with  the  Gaussian  radial  profile,  all  sub¬ 
sequent  simulations  were  initialized  with  the  Gaussian 
vortex,  as  described  in  Cl. 


d.  MM5  setup 

MM5  is  a  limited-area,  nonhydrostatic,  terrain-following, 
sigma-coordinate  numerical  model.  Grell  et  al.  (1994) 
provide  a  complete  description  of  MM5.  For  each  ex¬ 
periment,  a  coarse-resolution  simulation  was  first  per¬ 
formed  using  a  90-km  domain  with  a  size  of  approximately 
12  000  km  X  5000  km  in  the  zonal  and  meridional  di¬ 
rections,  respectively,  and  23  vertical  levels.  The  coarse- 
resolution  simulation  settings  included  full  Coriolis  force, 
the  simple  ice  (Dudhia)  moisture  scheme,  the  Grell  et  al. 
(1994)  cumulus  parameterization,  and  the  Medium-Range 
Forecast  (MRF)  model  boundary  layer  parameteriza¬ 
tion.  The  lateral  boundary  condition  for  all  runs  was  that 


Table  1.  Selected  parameters  for  basic  state  configuration. 


Shearcent 

CAPEeenl 

CAPEmax 

Gent  (K) 

(m  s~^  km^^) 

Widthez  (°) 

(j  kg-‘) 

(J  kg-') 

CA00SH15BZ60 

293 

1.5 

60 

0 

0 

CA22SH15BZ60 

303 

1.5 

60 

2247 

4200 

CA15SH30BZ60 

300 

3.0 

60 

1553 

4000 

CA15SH30BZ30 

300 

3.0 

30 

1553 

4100 

CA22SH08BZ60 

303 

0.75 

60 

2247 

4100 
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Table  2.  Parameters  of  inner  domains. 


Simulation 

Inner  domain  start  time  (h) 

Move  frequency  (h) 

Grid  dimensions* 

Grid  3 

Grid  4 

Grid  2 

Grid  3 

Grid  4 

CA00SH15BZ60 

120 

— 

24 

94  X  76 

229  X  172 

271  X  211 

CA22SH15BZ60 

120 

72 

8 

154  X  94 

289  X  199 

289  X  253 

CA15SH30BZ60-1 

42 

24 

8 

2II  X  109 

286  X  286 

286  X  286 

CA15SH30BZ60-2 

72 

24 

8 

211  X  109 

286  X  286 

286  X  286 

CA15SH30BZ30 

24 

24 

8 

211  X  109 

286  X  286 

286  X  286 

CA22SH08BZ60 

192 

— 

24 

211  X  109 

286  X  286 

286  X  286 

*  Grids  2,  3,  and  4  are  the  30-,  10-,  and  3.3-km  grids,  respectively. 


of  the  basic-state  environment.  As  in  Cl,  the  surface  fluxes 
of  latent  and  sensible  heat  as  well  as  radiative  processes 
were  neglected,  but  the  surface  flux  of  momentum  was 
retained. 

The  coarse-resolution  simulation  was  integrated  for¬ 
ward  in  time  to  provide  an  estimate  of  when  deep  con¬ 
vection  might  initiate.  This  estimate  was  then  used  for 
purposes  of  guiding  the  placement  and  movement  of 
finer-resolution  grids  for  a  second,  higher-resolution 
simulation.  To  provide  for  a  near-cloud-resolving  simu¬ 
lation  of  the  convective  processes  associated  with  the 
MCV  growth,  three  inner  nested  domains  were  used  with 
grid  intervals  of  30, 10,  and  3.3-km.  While  the  finest  grid 
interval  of  3.3  km  is  not  expected  to  fully  represent  pro¬ 
cesses  on  the  cloud  scale  (Bryan  et  al.  2003),  we  expect  to 
resolve  those  processes  that  are  essential  to  the  genera¬ 
tion  of  MCVs.  The  Grell  et  al.  (1994)  cumulus  parame¬ 
terization  was  used  on  the  90-,  30-,  and  10-km  domains 
but  turned  off  in  the  3.3-km  grid.  The  MRF  PBL  scheme 
was  used  on  the  30-km  domain,  and  the  Blackadar  PBL 
scheme  was  used  on  the  10-  and  3.3-km  domains. 

In  this  finer-resolution  simulation,  the  inner  grids  were 
initialized  when  a  contiguous  area  of  convective  accu¬ 
mulated  precipitation  developed  in  the  coarse-resolution 
simulation,  grew  over  at  least  two  consecutive  12-h  in¬ 
tervals  (e.g.,  12  h,  24  h,  36  h,  etc.),  and  reached  values  of 
at  least  5  mm.  The  inner  grids  were  then  turned  on  at  the 
beginning  of  the  nearest  6-h  period  prior  to  this  onset. 
For  example,  if  persistent,  convective  accumulated  pre¬ 
cipitation  was  first  observed  at  /  =  48  h  in  the  simulation, 
the  inner  grids  were  initialized  at  /  =  42  h  in  the  fine- 
resolution  simulation. 

Because  of  the  varying  conditions  of  shear  and  CAPE, 
it  was  necessary  to  move  the  innermost  two  domains  at 
different  intervals  for  each  experiment  as  indicated  in 
Table  2.  The  simulation  was  integrated  forward  in  time 
until  the  convective  system  grew  beyond  the  bounds  of 
the  3.3-km  domain. 

e.  Dry  simulations 

In  addition  to  the  experiments  above,  we  conducted 
several  dry  simulations  to  investigate  the  effect  of  the 


latent  heating,  provided  by  the  convective  development, 
on  the  moist  baroclinic  system.  In  these  simulations,  the 
corresponding  moist  case  was  run  long  enough  for  the 
deep  convection  to  initiate,  the  data  were  saved,  and 
then  the  runs  restarted  with  the  latent  heating  turned 
off.  For  the  CA22SH15BZ60  simulation,  a  dry  restart 
was  also  performed  several  hours  before  the  initiation  of 
deep,  moist  convection.  A  summary  of  the  dry  runs  is 
provided  in  Table  3. 

f.  Analysis  tools 

1)  Energy  conversion  calculations 

As  in  Cl,  we  calculated  the  diabatic  (G^)  and  baro¬ 
clinic  (Ca)  conversion  terms  to  compare  the  energetics 
of  the  CA00SH15BZ60  and  CA22SH15BZ60  simula¬ 
tions.  These  conversion  terms  represent  the  baroclinic 
conversion  of  mean-state  APE  and  the  diabatic  conver¬ 
sion  of  latent  heat  energy,  respectively,  to  eddy  APE.  The 
terms  are  defined  by  the  following  formulas  (Norquist 
et  al.  1977): 


Ge  = 


'P2Q'T' 

— ^  dp, 

n  C  a 
Pi  p 


(1) 


c. 


~P2  [v'T’]d[T] 


dp 


‘P2  [(o'T']d[T]* 

P,  dp 


dp. 


(2) 


where  [(•)]  represents  a  zonal  average  of  the  quantity  (■) 
and  [(■)]  represents  an  area  mean.  A  prime  indicates  the 
deviation  from  the  zonal  average,  and  an  asterisk  in¬ 
dicates  a  deviation  of  the  zonal  average  from  the  area 
mean;  Q  is  the  diabatic  heating  rate  (J  kg^^  s^'),  T  is 
temperature,  Cp  is  the  heat  capacity  at  constant  pressure 


Table  3.  Restart  times  for  the  dry  simulations. 


Experiment 

Restart 
time  (h) 

Time  relative  to  deep 
convection  initiation  (h) 

CA22SH15BZ60 

140 

-10 

CA22SH15BZ60 

160 

10 

(second  dry) 

CA15SH30BZ30 

48 

4 
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for  dry  air  (J  kg  ^  K  ^),  and  cr  is  the  area-average  static 
stability,  where 


a 


p  ^mV 

R  dp  y 


also,  V  is  the  meridional  component  of  velocity,  w  is  the 
vertical  velocity  in  pressure  coordinates,  and  p\  and  p^ 
are  the  pressure  limits  of  integration. 

We  conducted  a  limited  area  analysis  on  the  90-km 
domain,  centered  on  the  surface  pressure  minimum,  and 
extending  15°  of  longitude  (approximately  1350  km) 
east  and  west  of  the  minimum  and  10°  (approximately 
900  km)  north  and  south  of  the  minimum.  The  analysis 
was  performed  on  the  90-km  domain  in  order  to  en¬ 
compass  the  entire  large-scale  system. 

2)  Vertical  shear  analysis 

Using  the  method  outlined  in  Davis  and  Trier  (2007), 
we  calculated  the  area  average  of  1^  and  1-8  km  ver¬ 
tical  shear  within  300  km  of  the  MCV  center.  We  per¬ 
formed  the  calculations  on  the  30-km  domain  and 
defined  the  MCV  center  to  be  the  location  of  the  mini¬ 
mum  700-hPa  geopotential  height  perturbation  from  the 
zonal  average.  The  zonal  and  meridional  wind  compo¬ 
nents  u  and  V  at  1,  4,  and  8  km  AGL  were  horizontally 
averaged  over  an  area  within  300  km  radius  of  the 
center.  The  vector  differences  among  these  horizontal 
averages  were  then  calculated. 


3.  Results 

The  timeline  on  which  we  present  the  results  of  the 
simulations  begins  upon  the  initiation  of  deep  moist 
convection  and  lasts  for  40  h  after  initiation.  We  identify 
the  initiation  to  be  the  time  when  the  first  model-derived 
reflectivity  greater  than  5  dBZ  appears  at  the  700-hPa 
level,  followed  by  growth  to  become  an  MCS. 

a.  General  overview  of  simulations 

Initiation  of  deep  moist  convection  in  the 
CA22SH15BZ60  simulation  occurs  approximately  150  h 
into  the  simulation.  Two  hours  later  (Fig.  2a),  the  con¬ 
vection  has  grown  into  a  small  convective  cluster.  Six 
hours  after  initiation  (Fig.  2b),  a  rather  symmetric  MCS 
has  formed,  featuring  enhanced  700-hPa  rear  inflow 
along  52°N  and  a  leading  convective  line  spreading  in 
a  horseshoe  pattern  around  this  rear  inflow.  An  area  of 
stratiform  precipitation  is  developing  from  the  rear  in¬ 
flow  into  the  center  of  the  MCS.  It  is  in  Fig.  2b  that  the 
cyclonic  circulation  in  the  northern  portion  of  the  MCS 
starts  to  become  visible.  Four  hours  later  (Fig.  2c),  a  well- 
developed  MCV  is  present  in  the  stratiform  precipitation 


area  behind  the  leading  convective  line.  This  leading  line 
extends  southwestward  from  the  developing  circulation 
but  also  wraps  cyclonically  around  the  north  side  of  the 
circulation.  Eighteen  hours  after  initiation  (Fig.  2d),  the 
same  basic  structure  is  present,  but  it  has  grown  upscale  to 
the  point  that  the  convection  has  expanded  outside  of  the 
3.3-km  grid.  The  aerial  extent  of  the  700-hPa  circulation 
has  grown  beyond  the  limits  of  the  innermost  grid. 

A  look  at  larger  scales  (Fig.  3)  reveals  a  close  associa¬ 
tion  between  the  MCS  and  a  deepening  moist  baroclinic 
cyclone.  The  development  of  deep,  moist  convection 
occurs  just  to  the  east  of  the  surface  pressure  minimum 
(Fig.  3a),  which  is  a  location  favored  for  organized, 
persistent  lifting  motion  in  moist  baroclinic  cyclones 
(Bluestein  1992;  Holton  2004).  Ten  hours  after  initia¬ 
tion,  the  large-scale  cyclone  is  deepening,  and  the  MCS 
is  spreading  outward  around  the  center  of  the  cyclone 
and  along  the  leading  edge  of  its  expanding  cold  pool. 
The  lack  of  surface  fluxes  in  the  simulations  allows  the 
cold  pool  development  to  be  readily  seen.  For  the  next 
20  h  (Figs.  3c  and  3d),  the  surface  pressure  minimum 
continues  to  deepen  as  the  MCS  spreads  at  the  leading 
edge  of  its  cold  pool  and  also  wraps  cyclonically  around 
the  surface  low.  Thirty  hours  after  initiation  (Fig.  3d), 
a  substantial  portion  of  the  convection  is  located  to  the 
north  of  the  surface  low.  This  location  is  influenced  by 
the  presence  of  moderate  CAPE  in  the  initial  state  north 
of  the  center  of  the  baroclinic  zone,  and  it  may  be  con¬ 
sidered  rather  atypical  of  real  atmospheric  MCS  evolution. 
This  result  motivated  us  to  carry  out  additional  experi¬ 
ments  with  a  narrower  baroclinic  zone  and  less  CAPE  to 
the  north  of  the  center  of  the  baroclinic  zone.  The  nar¬ 
rower  baroclinic  zone  provides  conditions  more  typical 
of  the  warm  season  with  less  basic-state  APE. 

An  examination  of  the  2-m  temperature  and  surface 
pressure  evolution  prior  to  initiation  (not  shown)  re¬ 
veals  that  the  MCS  is  forced  not  by  the  warm  core  vortex 
but  by  the  subsequent,  larger-scale  baroclinic  wave  ex¬ 
cited  by  the  initial  vortex.  The  larger-scale  development 
has  a  wavelength  of  roughly  3000  km,  determined  by 
measuring  the  half-wave  distance  from  the  surface  pres¬ 
sure  minimum  to  the  downstream  surface  pressure  max¬ 
imum  (Fig.  3a).  Based  on  the  parameters  in  the  center  of 
the  baroclinic  zone  of  the  basic  state  of  the  simulation 
(see  Fig.  1),  we  calculated  the  Rossby  deformation  radius 
NHIf  (where  N  is  the  Brunt-Vaisala  frequency,  H  = 
6000  m  is  the  scale  height  of  the  atmosphere,  and  /  the 
Coriolis  parameter)  to  be  700  km.  The  scale  of  devel¬ 
opment  is  therefore  within  the  fast-growing,  smaller 
wavenumber  portion  of  the  large  cyclone  growth  regime 
presented  in  Fig.  1  of  MM04. 

The  importance  of  the  larger,  baroclinic  system  is 
underscored  by  the  fact  that  in  every  experiment,  the 
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Fig.  2.  Model-derived  reflectivity  (dBZ;  shading)  and  horizontal  flow  (wind  barbs;  full  barb  =  5  m  s~^)  at  700  hPa  on  the  3.3-km  grid  at 
the  following  times  relative  to  initiation  of  deep  convection  in  the  CA22SH15BZ60  simulation:  (a)  2,  (b)  6,  (c)  10,  and  (d)  18  h.  The  top 
and  right  axes  indicate  longitude  and  latitude,  respectively;  the  bottom  and  left  axes  indicate  grid  indices  on  the  3.3-km  domain.  The 
equivalent  radar  reflectivity  factor  for  rain  (mm^  m~'^)  is  computed  following  Fovell  and  Ogura  (1988)  as  Z^  =  12^akNf,Xr'’(plpy^)^, 
where  a  is  the  nondimensional  ratio  of  backscattering  for  the  reflecting  substance  (1  for  water,  0.213  for  ice),  k  =  10^^  is  the  conversion 
factor  from  m^  to  the  units  of  Z^,  No  is  the  intercept  parameter  of  the  raindrop  size  distribution,  p  the  density  of  the  substance  of  interest 
(water  or  ice;  p^  is  the  density  of  water),  and  A  is  the  slope  of  the  raindrop  size  distribution:  A  =  {'TTp^NolpaQia)^’'^^  where  p^  is  the  density  of 
dry  air  and  the  rainwater  mixing  ratio.  The  reflectivity  is  expressed  in  decibels,  or  dBZ,  where  dBZ  =  10  logio(Ze). 


MCS  development  occurs  just  to  the  east  of  a  minimum  in 
surface  pressure  (Fig.  4).  The  location,  in  every  circum¬ 
stance  of  initiation,  is  in  the  typically  favored  area  for  lift, 
east  of  the  surface  low,  where  the  large-scale  flow  induces 
relatively  strong  and  persistent  low-level  frontogenesis 
(Bluestein  1992;  Hoskins  1976).  In  the  CA00SH15BZ60 
simulation  (Fig.  4a),  initiation  occurs  at  the  north  end  of 
this  zone  at  a  location  that  has  CAPE  of  only  300  J  kg~^. 


In  this  experiment,  the  surface  low  development  is  col¬ 
located  with  the  initial  perturbation  vortex. 

In  the  CA22SH15BZ60  simulation  (Fig.  4b),  the 
frontogenesis  is  somewhat  stronger,  and  at  the  point  of 
initiation,  CAPE  is  2200  J  kg~\  Farther  east  of  the  sur¬ 
face  low,  around  11°W,  is  a  second  maximum  in  CAPE 
that  is  a  remnant  from  the  initial  perturbation  vortex.  No 
MCS  ever  forms  near  this  second  maximum.  A  second 
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Fig.  3.  Sea  level  pressure  (black  lines;  contour  interval  2  hPa),  2-m  temperature  (blue  dashed  lines;  contour  interval  2  K),  and  model- 
derived  reflectivity  (dBZ;  shading)  on  the  30-km  grid  at  selected  times,  relative  to  initiation  of  deep  convection,  during  the 
CA22SH15BZ60  simulation:  (a)  2,  (b)  10,  (c)  18,  and  (d)  30  h.  The  top  and  right  axes  indicate  longitude  and  latitude,  respectively;  the 
bottom  and  left  axes  indicate  grid  indices  on  the  30-km  domain. 


area  of  strong  frontogenesis  to  the  south  of  the  surface 
pressure  minimum  is  associated  with  an  area  of  convec¬ 
tion  that  does  not  develop  an  MCV. 

Two  separate  MCSs  developed  in  the  coarse-resolution 
CA15SH30BZ60  simulation,  so  we  conducted  two  separate 
finer-resolution  simulations  for  this  experiment.  We  des¬ 
ignated  these  two  simulations  as  CA15SH30BZ60-1  and 
CA15SH30BZ60-2.  The  first  MCS  (CA15SH30BZ60-1; 


Fig.  4c)  developed  near  the  center  of  the  baroclinic  zone, 
and  the  second  (CA15S1T30BZ60-2;  Fig.  4d)  developed 
20  h  later  and  farther  southwest. 

In  the  CA15SF130BZ60-1  simulation  (Fig.  4c),  the 
strongest  area  of  frontogenesis  near  42°N,  10°W  is  at  the 
southern  end  of  an  area  of  less  organized,  weaker  con¬ 
vection  associated  with  the  initial  perturbation  vortex. 
This  earlier  convection  weakened  as  the  baroclinic 


Fig.  4.  CAPE  (J  kg^^;  shading),  surface  pressure  (hPa;  thick  black  contour  lines;  contour  interval  1  hPa),  sum  of  the  deformation  and 
divergence  terms  of  the  Miller  frontogenesis  function  [K  (100  km  h)^*;  thin  black  contour  lines],  and  horizontal  flow  (wind  barbs;  full  barb  = 
5  m  s^')  on  the  90-km  (outermost)  grid  immediately  preceding  initiation  of  deep  convection  for  the  following  experiments:  (a) 
CA00SH15BZ60;  (b)  CA22SH15BZ60;  (c)  CA15SH30BZ60-1;  (d)  CA15SH30BZ60-2;  and  (e)  CA15SH30BZ30.  The  frontogenesis 
contours  have  the  following  properties  in  the  individual  panels  [zero  contour  omitted;  units  of  K  (100  km  h)^*]:  (a)  low  0.01,  high  0.04, 
interval  0.01;  (b)  low  —0.08,  high  0.09,  interval  0.01;  (c)  low  0.06,  high  0.24,  interval  0.02;  (d)  low  —0.3,  high  0.5,  interval  0.1;  and  (e)  low 
-0.02,  high  0.03,  interval  0.01.  The  “X”  marks  the  location  of  the  initiation  of  deep  convection.  The  top  and  right  axes  indicate  longitude 
and  latitude,  respectively;  the  bottom  and  left  axes  indicate  grid  indices  on  the  90-km  domain. 
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cyclone  developed  farther  west,  near  30°W.  The  target 
MCS  then  initiated  within  this  baroclinic  cyclone,  un¬ 
associated  with  the  initial  vortex.  Such  an  outcome,  on 
its  own,  provides  strong  evidence  of  the  interdependency 
between  the  moist  baroclinic  cyclone  and  the  MCS. 

The  CA15SH30BZ60-2  MCS  (Fig.  4d)  develops  within 
the  weakest  cyclone  of  any  of  the  experiments  (in  terms 
of  surface  pressure  perturbation),  but  a  distinct  trough  is 
still  present.  The  weaker  surface  pressure  minimum  is 
partly  due  to  the  latitude  (37°N)  but  may  also  be  a  reflec¬ 
tion  of  its  development  south  of  the  strongest  baroclinic 
zone.  The  initiation  here  may  be  more  CAPE-driven  than 
cyclone-driven,  but  a  moist  baroclinic  cyclone  develops 
nonetheless. 

In  the  experiment  CA15SH30BZ30  (Fig.  4e),  the  baro¬ 
clinic  cyclone  develops  in  association  with  the  initial  per¬ 
turbation  vortex.  The  MCS  initiates  in  the  same  favored 
location  slightly  downshear  from  the  surface  pressure 
minimum  in  an  area  with  strong  meridional  variation  in 
CAPE. 

In  the  CA22SH08BZ60  experiment  (not  shown), 
neither  deep  moist  convection  nor  a  strong  surface  low 
ever  developed.  The  surface  pressure  minimum  attainted 
during  the  simulation  was  only  2  hPa  deeper  than  the 
basic  state.  The  shear  of  0.75  m  s^^  km^'^  in  this  par¬ 
ticular  experiment  was  likely  too  weak  to  be  considered 
baroclinically  unstable  since  it  is  only  one-fourth  the  value 
of  the  Eady-type  basic  state.  Without  any  significant  forc¬ 
ing  for  upward  vertical  motion  due  to  synoptic-scale 
processes,  no  organized  MCS  ever  developed. 

In  the  following  discussion,  we  estimate  the  strength  of 
the  MCV  by  taking  the  maximum  difference  in  zonal 
winds  between  the  northern  side  and  the  southern  side  of 
the  cyclonic  circulation  in  Fig.  5.  The  MCSs  in  all  simu¬ 
lations  (except  for  the  experiment  CA22SH08BZ60) 
evolve  similarly,  yet  there  are  some  noteworthy  dif¬ 
ferences.  Ten  hours  after  initiation  (Figs.  2c  and  5),  the 
simulations  all  feature  an  MCS  with  a  cyclonic  circu¬ 
lation  (MCV)  at  700  hPa,  but  the  size  of  the  MCS  and 
strength  of  the  circulation  vary  from  case  to  case.  In  the 
CA00SFI15BZ60  simulation  (Fig.  5a),  only  weak  CAPE 
(less  than  700  J  kg^^)  is  available  to  the  MCS,  and  con¬ 
sequently  the  size  of  the  MCS  and  the  strength  of  the 
circulation  are  much  less  than  in  the  other  conducted 
experiments.  Despite  the  strength  and  width  of  the  baro¬ 
clinic  zone  being  identical  to  that  in  the  experiment 
CA22SF115BZ60  (see  Fig.  2c),  the  MCS  is  much  less  ex¬ 
tensive,  and  the  circulation  is  much  weaker. 

In  the  CA15SH30BZ60-1  simulation  (Fig.  5b),  the 
MCS  and  MCV  have  stronger  700-hPa  winds  than  in  the 
CA00SF115BZ60  simulation  (Fig.  5a).  Because  of  mod¬ 
erate  CAPE  and  a  doubly  strong  baroclinic  zone  (see 
Table  1),  the  increased  wind  speed  is  expected.  Com¬ 


pared  to  the  CA22SHI5BZ60  simulation  (Fig.  3c),  how¬ 
ever,  the  convection  on  the  north  side  of  the  MCS  is 
somewhat  weaker,  and  the  MCV  is  only  equally  as  strong. 
Despite  the  stronger  baroclinic  zone  relative  to  the 
CA22SHI5BZ60  simulation,  the  lack  of  strength  can  be 
explained  by  two  factors.  First,  there  is  less  CAPE  in  the 
center  of  the  baroclinic  zone  (see  Table  1),  and  the  CAPE 
diminishes  to  zero  north  of  the  baroclinic  zone.  Second, 
the  vortex  that  was  inserted  into  the  initial  state  reduced 
the  CAPE  to  the  east  of  the  developing  MCS.  The  ther¬ 
mal  and  moisture  advection  by  the  north-to-south  flow  on 
the  west  side  of  the  vortex  was  the  primary  cause  of  di¬ 
minished  CAPE.  Additionally,  the  disorganized  area  of 
deep,  moist  convection  that  developed  in  association  with 
the  initial  vortex  consumed  a  small  amount  of  CAPE. 

The  largest  MCS  in  all  the  simulations  was  in  the 
CA15SH30BZ60-2  experiment  (Fig.  5c).  Because  it  de¬ 
veloped  within  an  extensive  area  of  large  CAPE,  convec¬ 
tion  was  able  to  spread  north  and  south  of  the  developing 
MCV,  and  so  the  shape  of  the  MCS  most  strongly  re¬ 
sembles  that  of  the  CA22SFII5BZ60  simulation.  The 
strong  700- hPa  rear  inflow,  reaching  37  m  s^'^near35.5°N 
43.5°W,  demonstrates  the  impact  of  a  doubly  strong  baro¬ 
clinic  zone  relative  to  the  CA22SH15BZ60  simulation. 

As  in  the  CAI5SF130BZ60-I  simulation,  the  effects 
of  reduced  CAPE  to  the  north  of  the  baroclinic  zone 
(only  510  J  kg^^  at  50°N  compared  to  1513  J  kg^^  in 
the  CA22SH15BZ60  simulation)  are  evident  in  the 
CA15SH30BZ30  simulation  (Fig.  5d),  where  the  con¬ 
vection  on  the  northern  side  of  the  MCS  is  somewhat 
less  extensive  than  in  the  CA22SH15BZ60  (Fig.  2c) 
and  CA15SH30BZ60-2  (Fig.  5c)  simulations.  The 
CA15SFI30BZ30  simulation  (Fig.  5d)  has  the  strongest 
vortex  in  Fig.  5  with  a  maximum  difference  in  zonal 
winds  (from  the  north  side  to  the  south  side  of  the  vor¬ 
tex)  of  approximately  45  m  s^^ 

The  MCS  development  within  the  large-scale  baro¬ 
clinic  cyclone  looks  qualitatively  very  similar  in  all  ex¬ 
periments:  the  deep  moist  convection  is  located  within 
or  around  the  deepening  baroclinic  cyclone,  as  already 
revealed  in  Fig.  3.  The  dynamically  significant  differ¬ 
ences  among  the  experiments  are  most  easily  revealed  by 
plotting  the  time  series  of  the  surface  pressure  minimum 
(Fig.  6)  relative  to  initiation  time.  The  CA00SHI5BZ60 
cyclone  experiences  the  least  deepening  of  all  the 
cases,  which  is  not  surprising  given  the  modest  shear  and 
near-zero  CAPE  of  the  initial  state.  The  experiment 
CA22SHI5BZ60  has  the  same  basic  state  shear,  but  its 
cyclone  deepening  is  considerably  faster.  The  only  major 
difference  between  these  two  simulations  is  the  initial 
CAPE,  strongly  indicating  that  the  latent  heat  release 
within  an  MCS  has  a  major  effect  on  the  larger-scale, 
baroclinic  cyclone  development. 
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Fig.  5.  Model-derived  reflectivity  (dBZ;  shading)  and  horizontal  flow  (wind  barbs;  full  barb  =  5  m  s^')  at  700  hPa  on  the  3.3-km  grid 
approximately  10  h  after  initiation  of  deep  convection  for  the  following  experiments:  (a)  CA00SH15BZ60,  (b)  CA15SH30BZ60-1,  (c) 
CA15SH30BZ60-2,  and  (d)  CA15SH30BZ30.  The  top  and  right  axes  indicate  longitude  and  latitude,  respectively;  the  bottom  and  left 
axes  indicate  grid  indices  on  the  3.3-km  domain. 
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The  CA15SH30BZ60-2  simulation  deepens  more 
slowly  than  all  but  the  CA00SH15BZ60  experiment. 
The  cold  pool  in  this  case  is  strongest  (2-m  temperature 
perturbation  up  to  18  K),  thereby  increasing  the  surface 
pressure.  Additionally,  the  cyclone  develops  10°  south  of 
the  center  of  the  initial  baroclinic  zone,  which  places  it  in 
a  less  favorable  location  with  respect  to  the  large-scale 
dynamics. 

The  remaining  two  cases  experience  a  nearly  equally 
rapid  deepening  that  is  consistent  with  what  one  might 
expect  for  cyclones  developing  near  the  center  of  the 
baroclinic  zone  in  moderate  shear.  The  deepening  in  the 


CA15SH30BZ30  experiment  slows  about  25  h  after 
initiation  as  the  cyclone  moves  north  of  the  baroclinic 
zone  and  occludes.  In  the  CA15SH30BZ60-1  simula¬ 
tion,  the  baroclinic  zone  is  wider.  Thus,  the  cyclone  can 
continue  to  extract  energy  from  the  basic  state,  and  it 
continues  to  deepen. 

b.  Energy  conversion  as  related  to  scale  interaction 

The  results  of  the  energy  conversion  calculations 
show  dramatically  different  energy  conversion  in  the 
CA22SH15BZ60  and  CA00SH15BZ60  simulations 
(Fig.  7),  which  is  consistent  with  the  differences  presented 


1130 


MONTHLY  WEATHER  REVIEW 


Volume  138 


Cu 


1010-  ^ 


im 

I  1000-1 

o 
u 
0. 

"3 

J  990-1 

R 
01 
y) 


980- 


-20 


- - 

V.  \ 

V"x 

S:-.. 

« • 

- CA00SH15BZ60  \ 

. CA22SH15BZ60 

. CA15SH30BZ60-1 

- CA15SH30BZ60-2 

- CAI5SH3()BZ30 

“• - 1 - «- 


0 


“T" 

20 


~1 

40 


^  "  ^inil 


Fig.  6.  Time  series  (h)  of  minimum  surface  pressure  (hPa;  taken 
from  90-km  grid)  relative  to  the  time  of  initiation  of  deep  convection 
(finit)  for  the  CA00SH15BZ60  (solid),  CA22SH15BZ60  (dashed), 
CA15SH30BZ60-1  (dotted),  CA15SH30BZ60-2  (dotted-dashedX 
and  CA15SH30BZ30  (triple  dotted-dashed)  simulations. 


in  Fig.  6.  The  CA22SH15BZ60  simulation  (Fig.  7a)  has 
more  than  an  order  of  magnitude  greater  diabatic  pro¬ 
duction  of  APE  {Ge)  than  the  CA00SH15BZ60  run 
does  (Fig.  7b).  The  evolution  of  Ge,  as  well  as  the 
conversion  of  basic-state  APE  to  eddy  APE  (Ca)  for  the 
CA22SH15BZ60  simulation,  shows  some  rather  in¬ 
teresting  characteristics.  First,  the  diabatic  term  Ge 
undergoes  a  sudden  increase  as  the  convection  develops, 


but  the  baroclinic  conversion  Ca  diminishes.  We  suspect 
that  during  this  period  some  kinematic  reorganization  is 
occurring  because  as  the  MCS  develops,  its  cold  pool 
spreads  throughout  a  portion  of  the  warm  sector  east  of 
the  surface  low,  interfering  with  the  baroclinic  conver¬ 
sion  process  of  the  large-scale  system  (see  Fig.  3).  Mean¬ 
while,  as  will  be  seen  in  section  3e,  the  MCV  undergoes 
rapid  intensification  (between  950  and  450  hPa),  and 
discernable  cyclonic  rotation  develops  at  the  MCV  scale 
(Figs.  2a-c).  Thereafter,  the  baroclinic  conversion  term 
Ca  increases  again  and  becomes  dramatically  larger  than 
in  the  CA00SFI15BZ60  simulation  (Fig.  7b).  Rapid  deep¬ 
ening  of  the  surface  low  resumes  16  h  after  initiation. 

In  both  experiments,  once  convection  develops  and 
the  larger-scale  system  undergoes  a  few  hours  of  ad¬ 
justment,  overall  system-scale  deepening  (which  can  be 
considered  a  combination  of  the  MCV  scale  and  the 
synoptic-scale  deepening)  occurs  much  more  quickly. 
The  deepening  is  far  greater  for  the  CA22SFI15BZ60 
experiment  than  in  the  CA00SF115BZ60  simulation. 

As  has  been  discussed  in  sections  3a  and  3b  and  in  Cl, 
the  MCS  development  depends  on  the  larger,  moist 
baroclinic  cyclone,  which  provides  the  deep,  persis¬ 
tent  lift  and  low-level  inflow  of  moisture  to  initiate  the 
convection.  The  results  of  the  energy  conversion  calcu¬ 
lations  here — particularly  with  the  more  rapid  deepen¬ 
ing  in  the  CA22SH15BZ60  experiment  relative  to  the 
CA00SH15BZ60  experiment — suggest  that  the  relation¬ 
ship  also  acts  conversely.  That  is,  the  MCV  has  a  feed¬ 
back  on  the  developing  baroclinic  system,  causing  it  to 
deepen  more  quickly  than  it  would  in  the  absence  of  the 
strong  diabatic  heating. 


Fig.  7.  Time  series  of  the  minimum  surface  pressure  (hPa;  dotted  line),  conversion  of  basic-state  APE  to  eddy  APE  {Ca’,  kg  s^^;  solid 
line),  and  conversion  of  diabatic  heat  energy  to  eddy  APE  (G^;  kg  s^"*;  dashed  line),  computed  on  the  90-km  grid  for  a  an  area  with 
dimensions  of  30°  longitude  X  20°  latitude  centered  on  the  surface  pressure  minimum  for  the  (a)  CA22SH15BZ60  and  (b) 
CA00SH15BZ60  simulations. 
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Fig.  8.  Energy  conversion  comparison  between  moist  simulations  (black  lines)  and  simulations  with  latent  heating  shut  off  (gray  lines)  for 
the  (a)  CA22SH15BZ60  and  (b)  CA15SH30BZ30  experiments.  For  notation,  see  Fig.  7. 


c.  Dry  simulations 

Despite  the  initial  decrease  of  baroclinic  conversion 
Ca  as  the  MCS  develops,  the  net  effect  of  the  diabatic 
heating  on  Ca  is  positive.  The  results  of  all  three  dry 
restarts  (see  Table  3)  show  that  the  baroclinic  conver¬ 
sion  Ca  is  less  whenever  the  diabatic  heating  is  removed 
from  the  simulation  (Fig.  8).  In  the  CA22SH15BZ60 
restarts  (Fig.  8a),  the  magnitude  of  Ca  depends  on  the 
time  the  restart  is  made.  If  the  restart  is  performed  10  h 
before  initiation,  Ca  reaches  only  one-third  its  corre¬ 
sponding  value  in  the  moist  simulation.  If  the  restart  is 
performed  10  h  after  initiation,  Ca  remains  nearly  as 
large  as  in  the  moist  case  until  approximately  t  —  finit  = 
20  h.  Thereafter,  it  remains  approximately  constant  while 
Ca  of  the  corresponding  moist  experiment  continues  to 
increase.  In  both  dry  restarts,  the  deepening  of  the  surface 
low  is  less  rapid  than  in  the  moist  simulation,  and  in 
neither  dry  case  does  the  deepening  rate  appear  expo¬ 
nential.  The  lack  of  apparent  exponential  deepening  in 
the  dry  restarts  does  not  necessarily  indicate  baroclinic 
stability.  Rather,  it  may  be  due  to  the  fact  that  the  “dry” 
e-folding  time  scale  is  long  compared  to  the  “moist” 
e-folding  time  scale,  making  it  difficult,  within  the  time 
scales  examined  here,  to  extract  the  unequivocal  dry 
exponential  growth  rate. 

In  the  CA15SFI30BZ30  dry  restart  (Fig.  8b),  the 
baroclinic  conversion  term  Ca  increases  about  half  as 
fast  as  it  does  in  the  corresponding  moist  experiment.  In 
the  moist  case,  Ca  reaches  its  peak  47  h  after  initiation 
and  then  begins  to  decrease  as  the  system  moves  north  of 
its  baroclinic  support  and  occludes.  In  the  dry  case,  the 
conversion  of  mean-state  APE  to  eddy  APE  is  slower,  so 
the  occlusion  process  is  delayed  and  the  system  con¬ 


tinues  deepening,  although  not  as  fast  as  in  the  moist 
case.  By  t  —  =  76  h,  Ca  in  the  dry  case  has  caught  up 

to  its  value  in  the  moist  case,  and  the  surface  low  is  about 
as  deep.  Given  these  results,  there  may  be  enough  baro¬ 
clinic  instability  in  this  case  for  the  rapid  deepening  to 
occur  without  convection.  It  is  nonetheless  evident  that 
the  convection  significantly  accelerates  the  intensification 
process. 

The  results  of  all  restarts  suggest  that  the  MCS  has 
a  projection  onto  a  packet  of  unstable  moist  baroclinic 
modes  and  promotes  baroclinic  development.  From 
a  qualitative  point  of  view,  the  overall  length  of  the 
large-scale  wave  changes  little  during  the  simulation 
(see  Fig.  3),  which  suggests  a  preferred  wavenumber 
near  the  most  rapidly  growing  moist  baroclinic  mode 
(MM04).  During  the  first  10  h  after  initiation  in  the 
CA22SH15BZ60  experiment,  the  MCV  intensifies  on  its 
own  scale  (essentially,  the  scale  of  the  convective  sys¬ 
tem),  but  later,  once  reorganization  of  the  system  oc¬ 
curs,  deepening  proceeds  on  the  larger  scales  as  well. 
These  same  results  are  seen  in  the  CA15SH30BZ30 
experiment.  Our  overall  results  lend  support  to  the 
findings  of  Galarneau  et  al.  (2009),  who  likewise  docu¬ 
mented  the  upscale  growth  of  an  MCV  and  its  sub¬ 
sequent  interaction  with  a  baroclinic  zone  from  10  to  14 
June  2003,  during  BAMEX.  Their  analysis  demon¬ 
strated  the  influence  of  MCV  dynamics  on  guiding  the 
organization  of  a  larger-scale  baroclinic  wave. 

d.  PV  structure 

The  initial  PV  production  is  quite  strong,  localized, 
and  associated  with  the  initial  convective  elements  (Figs. 
2a,c  and  9a,b).  Thereafter,  as  the  cold  pool  grows  and 
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Fig.  9.  Potential  vorticity  (PVU;  shading),  model-derived  reflectivity  (dBZ;  black  contour  lines;  contour  interval  of  10  dBZ  from  10  to 
40  dBZ),  and  horizontal  flow  (wind  barbs;  full  barb  =  5  m  s^^)  at  700  hPa  on  the  30-km  grid  for  the  CA22SH15BZ60  experiment  at  the 
following  times  relative  to  the  initiation  of  deep  convection:  (a)  2,  (b)  10,  (c)  18,  and  (d)  30  h.  The  top  and  right  axes  indicate  longitude  and 
latitude,  respectively;  the  bottom  and  left  axes  indicate  grid  indices  on  the  30-km  domain.  The  straight  lines  in  the  figure  panels  indicate  the 
positions  of  the  vertical  cross  sections  appearing  in  Fig.  10. 


the  convection  moves  radially  outward  from  the  MCV 
(Figs.  9c, d),  the  strongest  PV  production  moves  outward 
as  well.  Viewed  on  the  3.3-km  grid  (Conzemius  and 
Montgomery  2009),  the  leading  line  of  convection  pro¬ 
duces  strong  couplets  of  positive  and  negative  PV,  with 
magnitudes  reaching  30  PV  units.  These  centers  fill,  merge, 
and  are  advected  rearward  from  the  convection  into  the 
stratiform  region,  closer  to  the  MCV  center,  where  PV 
production  continues.  Davis  and  Galarneau  (2009)  docu¬ 
mented  a  similar  process  for  two  BAMEX  MCVs.  On  the 
coarser  30-km  grid  (Fig.  9),  the  process  results  in  a  strong, 
positive  PV  anomaly  appearing  behind  the  leading  con¬ 
vective  line.  Upscale  growth  of  the  entire  system  occurs  as 


the  PV  associated  with  the  continuing  convection  becomes 
cyclonically  wrapped  around  the  original  positive  PV 
perturbation.  Although  this  line  of  newly  created  PV 
wraps  back  into  the  center  of  circulation  (because  of  dif¬ 
ferential  advection),  once  it  is  formed  there  is  little  or  no 
shearing  or  filamentation  of  this  new  PV  that  would 
indicate  a  Rossby  wave  “spinup”  process  (Carr  and 
Williams  1989;  Montgomery  and  Kallenbach  1997). 

Figure  10  shows  the  PV  vertical  cross  sections  corre¬ 
sponding  to  the  lines  shown  in  Fig.  9.  The  cross  sections 
reveal  the  typical  MCV  structure,  with  strong  positive 
PV  anomalies  throughout  the  troposphere  and  negative 
anomalies  found  downshear  in  the  upper  troposphere 
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Fig.  10.  Potential  vorticity  (PVU,  10^*"  K  kg^^;  shading),  potential  temperature  (white  lines;  contour  interval  2  K),  and  me¬ 

ridional  wind  (black  lines;  contour  interval  2  m  s^*)  on  the  30-km  grid  in  the  CA22SH15BZ60  simulation  for  the  following  cross  sections 
and  times  relative  to  the  initiation  of  deep  convection:  (a)  2  h  from  51.4°N,  27.5°W  to  51.4°N,  7.5°W;  (b)  10  h  from  52.2°N,  24°W  to  52.2°N, 
4°W;  (c)  18  h  from  52°N,  22.5°W  to  52°N,  2.5°W;  and  (d)  30  h  from  53°N,  17.5°W  to  53°N,  2.5°E.  See  Fig.  9  for  the  horizontal  axes  of  these 
cross  sections. 


and  stratosphere.  The  vortex  at  700  hPa  rapidly  in¬ 
tensifies  during  the  first  10  h  after  initiation  (Figs.  10a, b) 
and  then  grows  upscale  while  more  slowly  increasing  in 
intensity  (Figs.  10c, d).  Such  an  evolution  of  vertical  PV 
structure  is  consistent  with  the  outward  movement  of 
the  diabatic  PV  production  associated  with  the  leading 
convection.  At  later  stages  in  the  development  of  the 
mesoscale  vortex  (Figs.  10c,d),  most  of  the  increase  in 
meridional  winds  is  associated  with  the  amplification  of 
the  large-scale  baroclinic  wave.  By  30  h  after  initiation 
(see  Fig.  lOd),  the  diameter  of  the  vortex,  defined  in 
terms  of  maximum  meridional  wind,  is  approximately 
750  km,  which  is  somewhat  more  consistent  with  the  size 
of  a  typical  baroclinic  wave  than  that  of  an  MCV  (Davis 


and  Trier  2007).  The  MCV  in  this  experiment  undergoes 
significant  upscale  growth  as  the  convection  is  continu¬ 
ally  adding  to  the  PV  within  the  entire  system.  Such 
upscale  growth  has  also  been  observed  in  situations 
when  repeated  convection,  occurring  over  several  days, 
allows  for  reintensification  of  the  MCV  (Fritsch  et  al. 
1994;  Galarneau  et  al.  2009). 

In  summary,  as  the  MCS  grows  in  size  and  the  leading 
convective  elements  advance  outward  (Fig.  9),  the  PV 
production  expands  outward  as  well,  and  the  whole  PV 
anomaly  undergoes  an  upscale  growth.  In  fact,  by  the 
end  of  the  simulation,  it  becomes  difficult  to  clearly  dis¬ 
tinguish  the  MCV  from  the  deepening  baroclinic  system 
as  if  the  two  were  involved  in  a  seamless  transition. 
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Fig.  11.  Potential  vorticity  (PVU;  shading),  temperature  (K;  black  contour  lines)  and  horizontal  flow  (wind  barbs;  full  barb  equals 
5  m  s^')  at  700  hPa  on  the  30-km  grid  for  the  following  experiments  approximately  10  h  after  the  initiation  of  deep  convection:  (a) 
CA00SH15BZ60;  (b)  CA15SH30BZ60-1;  (c)  CA15SH30BZ60-2;  and  (d)  CA15SH30BZ30.  The  top  and  right  axes  indicate  longitude  and 
latitude,  respectively;  the  bottom  and  left  axes  indicate  grid  indices  on  the  30-km  domain.  The  straight  lines  in  the  figure  panels  indicate  the 
positions  of  the  vertical  cross  sections  appearing  in  Fig.  12. 


Intense  PV  production  occurs  at  the  mesoscale  (along  the 
main  squall  line)  and  becomes  entrained  into  the  synoptic- 
scale  PV  anomaly.  This  is  similar  to  the  process  Trier  et  al. 
(1997)  have  found  for  the  entrainment  of  PV  into  meso¬ 
scale  line-end  vortices. 

Essentially  all  the  idealized  cases  experience  this 
characteristic  evolution.  By  10  h  after  initiation  (Figs.  11 
and  12),  all  feature  a  positive  PV  anomaly  immediately 
behind  the  line  of  heaviest  reflectivity,  and  all  have  the 
same  characteristic  PV  vertical  structure  (Cl;  Fritsch 
et  al.  1994;  Raymond  and  Jiang  1990). 

Despite  these  commonalities,  a  few  dynamically  im¬ 
portant  differences  arise  among  the  experiments.  The 


CA00SH15BZ60  run  (Figs.  11a  and  12a),  with  its  rela¬ 
tively  small  CAPE  and  therefore  small  latent  heat  con¬ 
tent,  has  the  least  diabatic  PV  production.  Its  MCV  is  the 
weakest  among  all  the  cases,  reaching  roughly  5  m  s^^ 
tangential  velocity  (we  take  the  maximum  tangential  ve¬ 
locity  as  the  average  absolute  value  of  maximum  and 
minimum  meridional  winds).  The  MCV  is  also  suffi¬ 
ciently  sheared  that  its  tilt  becomes  visible  at  these  scales 
(note  that  the  aspect  ratio  of  these  plots  does  not  easily 
reveal  tilt).  The  CA00SH15BZ60  experiment  also  has 
the  lowest  height  of  maximum  tangential  winds  (which  we 
identify  by  the  average  height  of  the  minimum  and  max¬ 
imum  meridional  wind  component)  at  about  840  hPa. 
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Fig.  12.  Vertical  cross  section  of  the  potential  vorticity  (PVU,  or  10^*’  s~*  K  kg^*;  shading),  potential  temperature  (white  contour 

lines;  interval  2  K),  and  meridional  velocity  (black  lines;  contour  interval  2  m  s^*;  dashed  lines  for  negative  contours)  approximately  10  h 
after  the  initiation  of  deep  convection  for  (a)  CA00SH15BZ60  run  along  a  cross  section  from  49.4°N,  13°W  to  49.4°N,  7°E;  (b) 
CA15SH30BZ60-1  from  48.4°N,  28°W  to  48.4°N,  10°W;  (c)  CA15SH30BZ60-2  from  36.3°N,  49.5°W  to  36.3°N,  35°W;  and  (d) 
CA15SH30BZ30  from  47.4°N,  39°W  to  47.4°N,  21°W.  See  Fig.  11  for  the  horizontal  axes  of  these  cross  sections. 


In  the  CA22SH15BZ60  simulation  (Figs.  9b  and  10b), 
the  MCV  is  much  stronger  and  is  somewhat  taller  than  in 
the  CA00SFtl5BZ60  simulation.  It  develops  a  stronger 
cold  pool,  and  its  maximum  tangential  winds  (14  m  s^') 
occur  near  800  hPa,  higher  than  in  the  CA00SH15BZ60 
run. 

The  greater  APE  in  the  basic  state  of  the 
CA15SH30BZ60  experiment  appears  to  have  a  strength¬ 
ening  effect  on  the  MCVs.  The  vortices  in  both 
CA15SH30BZ60  simulations  are  as  strong  as  or  stronger 
than  that  in  the  CA22SH15BZ60  simulation.  Despite 
access  to  less  CAPE,  the  MCV  in  the  CA15SH30BZ60-1 
simulation  (Figs.  11b  and  12b)  has  tangential  winds  of 
15  m  s^\  The  CA15SH30BZ60-2  MCV  (Figs.  11c  and 


12c),  with  its  access  to  undiminished  CAPE  and  significant 
shear,  has  tangential  winds  with  values  up  to  20  m  s^\ 
The  CA15SH30BZ60-2  vortex  also  has  the  highest 
altitude  of  maximum  tangential  winds.  Because  the 
CA15SF130BZ60-2  MCV  is  at  the  lowest  latitude  of  all 
the  experiments  (36.5°N),  the  thermal  vertical  structure 
of  the  basic  state  (higher  tropopause  height,  larger  CAPE) 
allows  for  much  more  vigorous  convection,  and  the  max¬ 
imum  heating  rate  occurs  at  a  higher  altitude.  This  point 
will  be  further  elucidated  in  section  3e. 

The  CA15SH30BZ30  simulation  (Figs,  lid  and 
12d)  has  at  least  as  strong  a  vortex  as  in  the  two 
CA15SH30BZ60  simulations.  The  shear  and  CAPE 
in  the  background  state  are  the  same  as  in  the 
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CA15SH30BZ60  experiments,  but  the  baroclinic  zone  is 
only  half  as  wide,  resulting  in  greater  cyclonic  relative 
vorticity  on  the  poleward  side  of  the  baroclinic  zone 
where  the  MCV  develops. 

e.  Further  comparison  with  BAMEX  observations 

An  initial  comparison  with  BAMEX  observations  was 
made  in  Cl.  Davis  and  Galarneau  (2009)  have  made 
more  detailed  comparisons  between  simulations  and 
observations  for  two  BAMEX  lOPs.  Our  intent  here  is 
to  perform  a  cursory  investigation  into  the  similarities 
and  differences  between  the  simulated  MCVs  in  the 
present  study  and  the  observed  MCVs  in  BAMEX.  We 
will  identify  the  likely  causes  of  any  differences. 

The  two  most  important  differences  between  the 
simulations  and  the  BAMEX  lOPs  were  that  (i)  the 
simulations  were  performed  in  an  environment  with 
relatively  little  initial  baroclinic  wave  amplitude  (our 
initial  condition  is  a  sum  of  surface  edge  waves  in  the 
Eady  model)  compared  to  typical  atmospheric  cases 
and,  therefore,  nearly  pristine  mean-state  APE;  and  (ii) 
except  for  the  CA00SH15BZ60  experiment,  the  initial 
state  had  untapped  areas  of  CAPE  in  the  center  of  the 
baroclinic  zone  and  everywhere  equatorward.  Because 
of  the  extensive  areas  of  CAPE  and  largely  untapped 
APE  in  the  initial  state,  the  simulated  MCSs  can  expand 
to  cover  exceptionally  large  areas,  and  their  MCVs,  with 
ever-present  convection,  continually  grow  in  size  to  the 
full  cyclone  wave  scale  over  the  course  of  the  simulation, 
as  both  CAPE  and  basic  state  APE  are  consumed.  In 
actual  atmospheric  cases,  convective  activity  is  intermit¬ 
tent,  thus  frustrating  the  upscale  growth  of  the  MCV  or 
allowing  the  MCV  to  be  destroyed  by  the  background 
shear  (Cl;  MM05;  Fritsch  et  al.  1994). 

The  comparisons  with  BAMEX  lOPs  are  probably 
best  made  about  8  to  16  h  after  convective  initiation  in 
the  idealized  simulations,  when  the  MCVs  are  relatively 
young.  After  that  point,  the  continued  expansion  of  the 
convective  systems  and  deepening  of  the  baroclinic  waves, 
although  not  entirely  unrealistic,  seem  less  typical  of  sum¬ 
merlike  conditions.  Table  4  shows  the  comparisons. 

One  particular  case  of  interest  was  lOP  8.  Its  MCV 
was  the  strongest  observed  in  BAMEX,  yet  the  shear 
over  its  depth  was  the  weakest.  In  an  attempt  to  sim¬ 
ulate  an  MCV  in  similar  conditions,  we  created  the 
CA22SH08BZ60  experiment,  but  organized  convection 
never  developed  in  that  simulation.  Upper-air  analyses 
for  lOP  8  (Davis  and  Trier  2007;  Davis  and  Galarneau 
2009)  actually  showed  the  target  MCV  to  be  on  the  cy¬ 
clonic  shear  side  of  a  band  of  roughly  20  m  s^'  mid- 
tropospheric  flow.  The  area  just  south  of  the  MCV  had 
larger  vertical  shear  in  a  narrow  baroclinic  zone,  which  is 
more  typical  of  warm  season  conditions  than  the  wider 


baroclinic  zones  of  the  idealized  simulations  in  the  pres¬ 
ent  study. 

For  all  simulated  MCVs  in  Table  4,  the  vertical  shear 
is  larger  than  in  lOP  8.  Some  insight  into  the  reasons 
for  this  can  be  gained  by  plotting  the  time  series  of 
shear  (see  section  2f)  for  the  CA22SH15BZ60  and 
CA15SH30BZ30  simulations  (Fig.  13).  Although  the 
basic  state  shear  was  weaker  in  the  CA22SH15BZ60 
experiment  (Fig.  13a),  the  shear  increased  during  the 
simulation  (Fig.  13a).  A  closer  look  at  the  8-km  winds 
and  PV  structure  in  the  CA22SH15BZ60  simulation 
(not  shown)  revealed  that  the  upper  tropospheric  neg¬ 
ative  PV  anomaly  was  displaced  downshear  from  the 
midtropospheric  PV  positive  anomaly.  This  displace¬ 
ment  resulted  in  increased  southwesterly  flow  directly 
above  the  low-level  vortex  at  z  =  1  km.  In  essence,  the 
shear  remained  strong  because  the  MCV  remained  in 
the  baroclinic  region  throughout  its  lifetime. 

In  the  CA15SH30BZ30  experiment  (Fig.  13b),  the 
narrower  baroclinic  zone  allowed  the  growing  MCV  to 
move  north  into  a  region  with  less  shear.  Although  (for 
the  sake  of  consistency  in  our  analysis)  the  shear  values 
reported  in  Table  4  are  significantly  larger  than  the  values 
in  BAMEX  lOP  8,  the  shear  in  the  simulated  MCV 
eventually  decreased  to  values  consistent  with  lOP  8.  By 
that  time,  however,  the  system  scale  was  much  larger  than 
the  MCV  reported  in  lOP  8.  In  all  other  simulations,  the 
broad  baroclinic  zone  did  not  allow  the  MCV  to  move  out 
of  the  region  of  deep  layer  shear,  so  the  reported  values 
changed  less  than  in  the  CA15SH30BZ30  experiment. 

In  general,  the  shear  in  the  simulations  covers  the  low  to 
middle  range  of  the  BAMEX  lOPs.  Further  simulations 
with  relatively  more  intense  baroclinic  zones  could  be 
conducted  to  explore  the  upper  range  of  BAMEX  shear. 
Because  the  barochnic  zone  in  the  CA15SH30BZ30  ex¬ 
periment  was  still  somewhat  wider  than  that  prior  to 
lOP  8,  further  simulations  could  also  employ  yet  nar¬ 
rower  initial  baroclinic  zones,  but  it  is  already  clear  from 
the  present  simulations  that  if  the  MCV  is  able  to  move 
out  of  the  region  of  strong  shear,  it  is  more  likely  to 
acquire  characteristics  similar  to  the  lOP  8  MCV. 

Overall,  the  strength  of  the  MCVs  produced  in  the 
simulations  is  in  rather  good  agreement  with  the  in¬ 
tensity  of  the  MCVs  observed  in  BAMEX.  The  very 
small  CAPE  of  the  CA00SH15BZ60  experiment  results 
in  a  weaker  MCV,  but  the  rest  of  the  experiments  de¬ 
velop  MCVs  that  are  roughly  as  intense  as  the  observed 
MCVs  (as  measured  by  V„),  with  tangential  winds 
generally  within  the  range  of  that  observed  in  BAMEX. 
The  CA15SH30BZ60-2  MCV  was  a  bit  stronger  at 
20  m  s^\  but  its  CAPE  also  lies  outside  the  range  ob¬ 
served  in  BAMEX.  About  8-12  h  after  their  devel¬ 
opment,  the  radii  of  maximum  winds  in  all  the  simulated 


April  2010 


CONZEMIUS  AND  MONTGOMERY 


1137 


Table  4.  Properties  of  the  simulated  MCVs  and  MCVs  observed  in  BAMEX;  is  the  maximum  tangential  velocity,  R„,  is  the  radius  of 
maximum  tangential  velocity,  and  P„,  is  the  pressure  level  of  maximum  tangential  velocity. 


Simulation 
or  lOP 

Date 

(2003) 

V,„  (m  s 

Rn,  (km)* 

P,„  (hPa) 

CAPE 
(J  kg  “')** 

Shear,  1-4  km 
(m  s-‘) 

Shear,  1-8  km 
(m  s-‘) 

CA00SH15BZ60 

— 

5 

100  T 

850 

300 

5 

12 

CA22SH15BZ60 

— 

14 

100  t 

800 

2200 

4 

13 

CA15SH30BZ60-1 

— 

15 

100  t 

650 

1700 

9 

16 

CA15SH30BZ60-2 

— 

20 

150  t 

580 

4000 

10 

14 

CA15SH30BZ30 

— 

17 

150  t 

630 

2800 

10 

18 

1 

24  May 

13 

100 

480 

1000 

12 

14 

4 

2  Jun 

15 

150 

650 

1000 

14 

27 

5 

5  Jun 

12 

100 

660 

2000 

12 

32 

8 

11  Jun 

15 

150 

660 

2000 

2 

5 

15 

29  Jun 

10 

100 

540 

1700 

8 

13 

*  The  up  arrow  in  the  R,„  column  is  an  indicator  that  is  increasing  with  time. 

**  CAPE  in  the  BAMEX  lOPs  has  been  calculated  from  the  Eta  model  analysis  of  CAPE  at  the  point  where  the  MCS  was  first  defined.  In 
the  simulations,  the  CAPE  value  is  taken  from  the  location  of  deep  convection  initiation. 


vortices  fall  within  the  range  observed  in  the  lOPs.  This 
result  provides  strong  indication  that  the  structures  be¬ 
ing  simulated  herein  are  useful  analogs  of  the  observed 
MCVs.  As  already  noted,  however,  because  of  the  ex¬ 
perimental  setup  the  vortices  in  the  idealized  simulations 
continue  to  grow  in  size  after  the  period  of  comparison. 

One  other  final  point  of  comparison  is  the  altitude 
of  maximum  tangential  winds.  In  BAMEX,  the  range 
presented  is  at  a  somewhat  higher  altitude  (lower  pres¬ 
sure)  than  in  the  simulations.  In  the  simulations,  the 
altitude  of  maximum  tangential  winds  appeared  to  be 
related  to  latitude  and  CAPE,  both  of  which  also  have 
an  impact  on  the  level  of  maximum  diabatic  heating. 
The  initial-state  thermal  profiles  at  lower  latitudes  are 
more  conducive  to  a  heating  maximum  at  a  higher  level. 
The  equation  for  the  material  rate  of  change  of  PV 
(Holton  2004)  implies  cyclonic  vorticity  production  below 


the  level  of  maximum  heating  (the  PV  production  to 
which  we  refer  is  for  parcels  in  a  Lagrangian  setting).  By 
the  invertibility  principle,  in  the  approximately  moist 
neutral  conditions  within  the  MCS,  most  of  the  cyclonic 
PV  anomaly  will  be  manifest  as  a  cyclonic  wind  anom¬ 
aly.  Thus,  as  the  level  of  maximum  heating  is  raised,  the 
level  of  maximum  tangential  wind  will  also  rise. 

Indeed,  the  CAI5SH30BZ60-2  experiment  has  the 
highest  level  of  maximum  tangential  winds.  All  other 
simulated  MCVs  occurred  north  of  45°N  where  the 
tropopause  height  and  CAPE  were  lower.  The  BAMEX 
MCVs  were  all  observed  south  of  45°N. 

4.  Summary  and  conclusions 

The  summary  and  conclusions  of  this  paper  can  be 
stated  in  terms  of  the  initial  goals  set  forth  in  section  1. 


f-4,„(hours)  f (hours) 


Fig.  13.  Time  dependence  of  the  vector  difference  in  horizontal  velocity  between  1  and  4  km 
MSL  (solid  lines)  and  between  1  and  8  km  MSL  (dashed  lines)  for  the  (a)  CA22SH15BZ60  and 
(b)  CA15SH30BZ30  simulations.  The  time-dependent  plots  are  taken  from  the  simulations, 
and  the  time-constant  lines  represent  the  shear  in  the  background  state. 
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1)  The  conducted  simulations  reveal  an  important  syn¬ 
ergistic  (two-way)  interaction  taking  place  between 
the  mesoscale  (MCS-MCV)  development  and  the 
synoptic-scale  baroclinic  waves.  The  results  lend 
strong  support  to  the  hypothesis  that  the  MCS  and 
subsequent  MCV  development  are  dependent  on  the 
moisture  flux  and  lift  provided  by  a  larger-scale  sys¬ 
tem.  In  the  simulations,  the  large-scale  system  was 
a  growing  baroclinic  wave,  but  in  the  warm  season 
atmosphere,  it  is  quite  plausible  that  more  station¬ 
ary  large-scale  structures  such  as  the  low-level  jet  on 
the  North  American  plains  may  also  play  a  role  in 
providing  this  moisture  flux,  where  such  structures 
intersect  baroclinic  regions.  In  the  opposite  sense, 
the  growth  rate  of  the  synoptic-scale  system  is  en¬ 
hanced  by  the  diabatic  production  of  PV  and  eddy 
available  potential  energy  (APE)  provided  by  the 
MCS.  Together,  these  findings  point  to  a  synergistic 
interaction  between  the  MCV  and  the  larger  synoptic- 
scale  flow. 

2)  The  spiral-shaped  PV  structures  appearing  in  the  x-y 
cross  sections  of  the  simulations  are  caused  by  in¬ 
tense  PV  production  occurring  in  and  just  behind  the 
leading  convective  line  (at  the  edge  of  the  cold  pool) 
that  advances  outward  from  the  MCV  center.  These 
new  areas  of  PV  then  rotate  cyclonically  around  the 
original  center  and  gradually  become  entrained  into, 
and  more  axisymmetric  with,  the  main  PV  center  as 
the  whole  system  grows  upscale. 

3)  Comparisons  with  the  MCVs  observed  in  the  Bow 
Echo  and  Mesoscale  Convective  Vortex  Experiment 
(BAMEX)  show  that  the  latitudinally  broad  baro¬ 
clinic  zones  in  the  constructed  basic  states  limited  the 
variability  of  vertical  shear  experienced  by  the  sim¬ 
ulated  MCVs  relative  to  the  MCVs  of  BAMEX. 
Also,  the  deepening  of  the  surface  lows  is  enhanced 
relative  to  what  might  be  expected  in  warm  season 
conditions.  This  is  the  result  of  the  pristine  initial 
state,  in  which  baroclinic  wave  amplitude  was  ini¬ 
tially  very  small.  The  basic  state  exhibited  instability 
to  moist  baroclinic  waves,  and  very  large  areas  of 
CAPE  were  available  to  the  first  system  that  devel¬ 
oped  in  the  simulation.  Nevertheless,  the  structures 
that  developed  had  features  similar  to  those  observed 
in  BAMEX,  with  the  radii  of  maximum  winds  and  the 
values  of  tangential  winds  being  within  the  range 
observed  in  BAMEX. 

In  the  light  of  making  further  comparisons  with 
BAMEX  MCVs,  additional  realism  may  be  provided  in 
future  simulations  by  allowing  for  narrower  baroclinic 
zones  than  have  been  simulated  in  the  present  study. 
This  would  permit  MCVs  to  move  in  and  out  of  regions 


of  deep  layer  shear  and  thereby  allow  us  to  explore  the 
effects  of  changing  shear  on  MCV  dynamics.  Finally,  the 
sensitivity  of  the  MCV  dynamics  to  the  diurnal  cycle  has 
not  been  actively  explored  in  these  simulations.  These 
complex  and  physically  relevant  issues  await  further 
investigation. 
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